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Abstract 

A 2-orbital t-J model over the square lattice that describes low-energy electronic excitations in 
iron-pnictide high-T c superconductors is analyzed with Schwinger-boson-slave-fermion mean field 
theory and by exact numerical diagonalization on a finite system. When inter-orbital hole hopping 
is suppressed, a quantum critical point (QCP) is identified that separates a commensurate spin- 
density wave (cSDW) state at strong Hund's rule coupling from a hidden half-metal state at weak 
Hund's rule coupling. Low-energy spinwaves that disperse anisotropically from cSDW momenta 
are predicted at the QCP. Nested Fermi surfaces similar to those observed experimentally in iron- 
pnictide materials are also predicted in such case. 
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I. INTRODUCTION 



High-temperature superconductivity in iron-pnictide materials is achieved by injecting 
charge carriers into stoichiometric parent compounds that show commensurate spin-density 



these _ have predo—y iron 3rf-orbitaJ characterBwH, ~* «* at two- 
dimensional (2D) hole-type and electron-type Fermi surface pockets centered, respectively, 
at zero 2D momentum and at the 2D cSDW momenta (/i/2a)x(y) {(| [7] js]. Here, a is the iron 
2D lattice constant. Density functional theory (DFT) calculations correctly account for the 
Fermi surfaces that are also observed in the cSDW[9], but they predict an ordered moment of 
approximately 2 Bohr magnetons (/x.b) that is large compared to measured values that can be 
as low as 0.3 /xb [2]. Frustrated 2D Heisenberg models that assume local magnetic moments 
at the iron atoms can account for the weak cSDW observed in parent compounds, on the 



wave (cSDW) order over a tetragonal lattice of iron atoms llyl The charge carriers in 



mr 



other hand 



10|(ll|(l2| 



HQ. 



They cannot predict the observed Fermi surface pockets, 



however, because local-moment magnets are Mott insulators. 

Below, we introduce mobile holes into a two-orbital local-moment model for frustrated 
m a g net 1Sm that success^ describe the weah eSDW nature of ._t.de SyS te ms Q. 
Exact diagonalization of one hole that hops over a square lattice of iron atoms in a single 
layer obtains a robust cSDW groundstate at 1/2 the cSDW momenta when Hund's rule is 
obeyed (cf. ref. [151]) . Both the exact and a mean-field analysis find that two Fermi surface 
hole pockets centered at zero 2D momentum can emerge when Hund's rule is violated, at 
sufficiently weak inter-orbital hopping of electrons. The ground state in question is the 
half-metal state with antiferromagnetic order across the 3dr x +i y \ z and 3dr x -i y \ z orbitals of 
the iron atom JjJ. Proximity to a quantum critical point (QCP) that separates the hidden 
half-metal at weak Hund's coupling from a cSDW metal 16[ at strong Hund's coupling 
results in weak cSDW order and in low-energy spin-wave excitations at cSDW wave numbers. 
The latter disperse in a manner that is consistent with inelastic neutron scattering (INS) 
measurements in iron-pnictide metals 13] 18j • The critical spinwaves in turn result in nested 
Fermi surface pockets that are centered at cSDW wave numbers, in agreement with angle- 
resolved photoemission (ARPES) on iron-pnictide metalsQ. These have mixed electron- 
type and hole-type character, however, which is a feature that has also been seen in ARPES 
studies of iron-pnictide metals js]. 
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II. 2-ORBITAL t-J MODEL: MEAN FIELD THEORY AND EXACT RESULTS 



We shall now show that a hidden half-metal groundstate emerges from the following two- 
orbital t-J model over the square lattice, where double occupancy at a site-orbital is strictly 
forbidden: 

_^ l r i 2 — ^ 

H = - Y Y H( t ?' /3 cl,a,^i,/3, S +h.C.) + - Jo Y Y S *,« + Y Y J? ,P Si,aSj,P+ Y Y J2^i, a -Sj,P- 

(1) 

Above, Si ta is the spin operator that acts on the spin so = 1/2 state of the 3d(x+i y ) z (d+) 
or 3d( x -iy) z (d—) orbital, a = or 1, in the iron atom at site i. The latter runs over the 
square lattice of iron atoms that make up an isolated layer. The former basis of d± orbitals 



is the least localized one (cf. ref. 



191]). which maximizes the Hund's coupling, — J . It 



therefore minimizes the Hund's rule exchange energy in the 2-orbital t-J model (pQ). Nearest 
neighbor and next-nearest neighbor Heisenberg exchange across the links and 
is controlled by the exchange coupling constants J" and J^ ,/3 ', respectively. These are 
necessarily isotropic over the d± orbital basis. Correlated hopping of an electron in orbital 
a to a neighboring unoccupied orbital (3 is controlled by the hopping matrix element The 
Heisenberg model that corresponds to (0Q) in the absence of charge carriers possesses a QCP 
at large Sq that separates a cSDW at strong Hund's coupling from a hidden ferromagnet that 
shows spin order at weak Hund's coupling if off-diagonal frustration exists 14{]: e.g. 

j] =0, J i > 0, and jj = J 2 > 0. Here the superscripts || and _L refer to the relationship 
between the orbital indices a and /3. Recent DFT calculations show that the superposition of 
direct ferromagnetic exchange with super-exchange across nearest-neighbor iron atoms can 
result in the cancellation jf = 4]. The remaining positive exchange coupling constants 
are assumed to be due to the super-exchange mechanism [10f] [20] . 

Let us now turn off inter-orbital hopping: t^ = 0. Notice that antiferromagnetic order 
across the d+ and d— orbitals then remains intact in the presence of mobile holes 14]. 
This is a classical picture for a half-metal groundstate that shows \d+\d- spin order. 
To describe it, we adopt the Schwinger-boson (b) slave-fermion (/) representation for the 

h,a,sfta and S i,a = (V 2 ) E s , s / fi,ab{ a>s Cr s , a >b i>a ,s> with 

(2) 



correlated electron 
the constraint 
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2l|: ~d 



2SQ — &i, a ,-|-&i,a,t + h,a,i.Ka,l + fi,afi,a 
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imposed at each site-orbital to exclude double occupancy. From here on we set % = 1. 



Following Arovas and Auerbach 22J , we next rotate about the spin y axis by an angle 7r 
on one of the antiferromagnetic sublattices in order to decouple the spins: b i ^ — > — ^2,0,4. 
and bifl^ — > &i,o,t- Again following these authors, we now define mean fields that are set 
by the pattern of ferromagnetic (||) versus antiferromagnetic (_L) pairs of neighboring spins: 
Qo = (bi, a , s bi,p,s) , QI(2) = { h \,a,s h i,a,s) and Qi {2) = (bi^sbj^g) for (next) nearest-neighbor 
links, where a 7^ (3. We add to that list the mean field P} = \{fl a fj,a) for nearest-neighbor 
hopping of holes within the same orbital. The mean-field approximation for the t-J model 
Hamiltonian ([T]) then has the form H MF = H [Q, P] + H b + Hj, where H [Q, P] consolidates 
the bilinear terms among the mean fields, where 

H b = 5 ££{n»(*)K(*)&.(*) + b s (-k)bl(-k)) + n ± (k)[bt(k)bi(-k) + b s (-k)b s (k))} 

Z k s 



is the Hamiltonian for free Schwinger bosons, with 



n,|(A0 = *A+ £ ^J^-4(j; ll gf + 24p|)[l-7i(k)]-4^ ll Ql[l- 72 (k)] 

n=0,l,2 

Q±(k) = _ c <* y: z n j' n x Qi ln (k), 

n=0,l,2 

and where Hf = J2k £ f{k)fKk)f(k) ^ s the Hamiltonian for free slave fermions, with 
£ f(k) = ^tiQilife) — Above, k = (ko, k) is the 3-momentum for these excitations, 
with corresponding destruction operators b s {k) = N~ l l 2 Y} a =o S« e^ fc ° a+k ' r ^6j >ai < j and f(k) = 
_1 ' 2 Ei=o E; e ! ^ (ia+k ' ri '/j iQ . Here ko = 0,tt represent even and odd superpositions of the 
(i± orbitals, while = 2Np e denotes the number of sites-orbitals on the square lattice of iron 
atoms. Also above, z = 1 and zi( 2 ) = 4, 7o(k) = 1 and 71(2) (k) = |(cos k x (+)a + cosky^a), 
where k± = k x ± k y , while 5X is the boson chemical potential that enforces the constraint 
against double occupancy (T5J) on average over the bulk of the system. The concentration of 
mobile holes per site-orbital, x, sets the chemical potential of the slave fermions, fi. Last, 
the effect of mobile holes on the Heisenberg spin-exchange is accounted for by the effective 



exchange coupling constants[21] J' = (1 — x) 2 J. 

The solution to the above mean field theory is achieved by making the standard Bo- 
goliubov transformation of the boson field, b s (k) = (cosh 6k)P s {k) + (sinh 8k)Pl{— k), with 
cosh26' = Q\\/tUb and sinh 29 = —tt±/tUb, where Ub = — fi^) 1 / 2 is the energy eigenvalue 
of the (/3) boson. Enforcing the constraint against double occuppancy (T5]) on average then 
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yields the principal mean field equation 22j 



s + \ - \x = Af- 1 '£(cosh26 k )(n B [cj b (k)} + ~), (3) 

" ~ k " 

where n B denotes the Bose-Einstein distribution. Ideal Bose-Einstein condensation (BEC) 
into the two lowest-energy states at k = occurs as temperature T — > 0, in which case 
5X 0. The remaining self-consistent equations for the Schwinger-boson mean fields are 

Ql = Ar- 1 ^ 7n (k)(cosh2^)(n B [o; fe (A;)] + i) and 
k 1 

k 1 

After comparison with ([3]), ideal BEC as T — > implies the unique value Q = Sq for all five 
of these mean fields in the large-so limit. Last, the self-consistent mean field equation for 
intra-orbital hole hopping is Pf = A/" -1 J2k 7i(k)^F^/(k)], where rip is the Fermi distribution 
function. We henceforth assume a hole band at low doping, t\ < and i C 1, which implies 
two degenerate circular Fermi surfaces centered at zero 2D momentum with Fermi wave 
vector kpa = (Aiix) 1 ! 2 . This yields the amplitude Pf = x/2 for intra-orbital hole hopping. 
Inspection of the spectrum for Schwinger bosons, Wf,(k), yields a spin gap at cSDW wave 
numbers (7r/a)x(y) equal to 

A cS dw = (1 - x) 2 (2s )[(4J 2 ± - J 0c )(Jo - Joc)] 1/2 , (4) 

where — J 0c = 2(Jf L — jf ) — 4j| — (1 — x)~ 2 Sq l 2t\x is the critical Hund's coupling at which 
A 

cSdw ~^ 0. Notice that intra-orbital hole hopping stabilizes the hidden half-metal state. 
We therefore propose (i) that the normal state of iron-pnictide superconductors is described 
by the present hidden half- metal state, and (ii) that the cSDW/superconductor transition 
that these systems commonly exhibit 2J 7| is controlled by the QCP at Hund's coupling — J 0c . 
The linear increase of — Jq c with the concentration of holes x implies a charge-carrier-poor 
cSDW and a charge-carrier- rich superconductor, which is consistent with experiment. 

Transverse dynamical spin correlations are obtained directly from the above Schwinger- 
boson-slave-fermion mean field theory. In particular, we have 

(S X (y)S x (y))\ kojk!U} = "(I - x) 2 [G b *G* b + {~)Fb * •F&*]U(0)+fc 0l k,U;, 

where iGb(k,co) = (b s (k,u)bl(k,u)) and iFb(k,co) = (b s (k,co)b s (—k, -co)) are the regular 
and the anomalous Greens functions for the Schwinger bosons, and where the notation 



f * g denotes a convolution in frequency and momentum. This yields an Auerbach-Arovas 
expression for the dynamical spin correlator at T > . It is easily evaluated in the zero- 
temperature limit, where ideal BEC of the Schwinger bosons into the doubly degenerate 
k = ground state occurs. It contributes to half of the net transverse spin correlator, which 
in this limit and at large So reads 

i(S± ■ S ± )| fclW = (1 - x) 2 s (Q + /n^ 2 ([u b (k) - u]- 1 + [cj b (k) + u}- 1 ). (5) 

Here, Q± = Q\\±Q±. The above dynamical spin correlator coincides with the transverse spin 
susceptibility, x±(k,u), in the present zero-temperature limit by the fluctuation-dissipation 
theorem. Observe now that fi_(7r,k) and O+(0, k) both vanish at k = in general, whi 



fi_(0, k) and fi+(vr, k) both vanish at cSDW wave vectors k = (7r/a)x(y) at the QCP [14]. 
The identity cot = (fi_fi + ) 1 ^ 2 then ultimately yields that spinwaves at zero 2D momen- 
tum disperse isotropically as w&(k) = t>o|k|, while those at cSDW wave numbers disperse 
anisotropically as 

u, 6 (k) = [vKh - n/a) 2 + v 2 (k t / lcSDW ) 2 + A 2 cSDW ] 1/2 . (6) 

(See fig. Here, k\ and kt are the longitudinal and the transverse components of k with 
respect to the cSDW wave number. The longitudinal spin-wave velocity is 

v = 2s a(l - x) 2 ([J 1 ± - 4(x) + 2(4 - j|)] • [~ J + 2J± + 2Ji}) 1 / 2 (7) 

and the anisotropy parameter is 

IcSdw = ([2(4' + Jt) + 4(x) + Jft/Wl + 4) - 4(x) - Jt}) 1 ' 2 , (8) 

which is greater than unity. Here jf (a?) = j| + (1 — x)~ 2 Sq l t\x. Study of the spectral weights 
in expression (jSJ) for oj) then yields that the former spinwaves at zero 2D momentum 

are hidden (fc = n), while that the latter spinwaves at cSDW momenta are observable 
(ko = 0). Figure [1^ displays x_i_(fc, w) at the QCP in the observable channel, k = 0, 
assuming off-diagonal magnetic frustration and a low concentration of mobile holes: jf =0, 



JjL > 0, 4 1 = 0.37^ = J 2 X , tj| = -SJj 1 , i| = and i = 0.01. Setting SoJ^ ~ 70 meV 
yields a successful fitfljj to spin- wave spectra obtained from INS on the superconductor 



BaFe2-a;Co :r As2 



181 ] . The critical spin- wave dispersion in fig. [T^i also notably shows a 



local maximum at the Neel momentum (ir/a, 7r/a), which agrees with INS on the parent 



compound CaFe2As2 [17J. Proximity to the QCP also naturally accounts for the low values 
of the magnetic moment associated with cSDW order (/j, c sdw) that are seen in iron-pnictide 
parent compounds by neutron diffraction [2]. Indeed, fig. [Tb displays exact results for one 
hole roaming over a 4 x 4 square lattice of iron atoms at the QCP, where ^ c sdw is a 
fraction of the maximum possible ordered moment achieved in the true ferromagnetic state, 
I^Fe = (33/31) 1 / 2 2 / ub. (See end of section.) Last, Eq. (jSJ) coincides with the large-so result 
for x±(k, oj) obtained by one of the authors in the hidden ferromagnetic state at x = jl4]. 

The electronic structure of the hidden half-metal state can also be obtained directly 
from the above Schwinger-boson-slave-fermion mean field theory. In particular, the electron 
propagator is given by the convolution of the propagator for Schwinger bosons with the 
propagator for slave fermions: iG(k,u>) = G b * G*f\k )UJ , where iGf(k,uj) = (f(k,co)p(k,uj)). 
A standard summation of Matsubara frequencies yields the expression 



<? L 



(cosh 9 q f ngMg) l + {<1 " f )] +(sinh g^"*^ + M ~ £/(g 



lo - u b (q) + £f(q - k) ' ' u + Lo b (q) + e f (q - k) 

(9) 

Ideal BEC of the Schwinger bosons at 2D momentum q = results in the following co- 
herent contribution to the electronic spectral function at zero temperature and at large Sq: 
lmG co h(k,u) = s 7i5[uj + £/(k)]. It reveals the two degenerate hole bands expected from 
the classical picture of a half-metal state with spin order and with no inter-orbital 

hopping[14]. The fermion contribution to lmG(k,u>) above represents incoherent excita- 
tions. They show a gap A c sdw (@J at cSDW momenta. Those originating from the second 
term in ([9]) are combinations of a hole with a spinwave, with a total energy that lies below 
the Fermi level. The incoherent contribution originating from the first term in (Q is the 
time-reversed counterpart, and it doesn't contribute at energies below the Fermi level in 
the zero-temperature limit. Last, the ratio of the incoherent spectral function integrated 
over momentum in the vicinity k = or (7r/a)x(y) compared to the coherent counterpart is 
Y,' k lraG inc (k,uj)/Y,k^Gcoh(k,uj) = (7/4tts ) • [(-u)n {l /(vo/a) 2 } as A cSDW ->■ 0, where -u 
measures how far in energy the hole lies below the Fermi level, and where 7 is the anisotropy 
parameter of the spinwave dispersion in question (see fig. [T^l) . 

We will now evaluate the former incoherent contribution to the spectral function in the 
large-so limit as A c sdw ~^ 0, at energies just below the Fermi level. The previous long- 
wavelength approximations for the spinwave dispersion near zero 2D momentum and near 
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cSDW momenta ([6]) are then valid. Also valid is the longwavelength approximation for the 
dispersion of the slave fermions, £/(k) = (2so)( — £i)(|k| 2 — k 2 F )a 2 . The imaginary part of 
the pole in the second term of Eq. enforces energy conservation, — oj = w&(q) + £/(q — 
k). The inequality Vp > Vq implies diffuse electron bands centered at momenta zero and 
(7r/a)x(y) that close out at an energy hp ■Vo/'y below the Fermi level. Figure [2] shows the net 
spectral function in these regions at the QCP. The incoherent contribution was evaluated 
in the thermodynamic limit by integrating the 5-function in energy over radial momentum 
analytically, and by performing the remaining angular integral numerically. Figure^ shows 
that the electronic structure at zero 2D momentum is predominantly hole-type because of 
the coherent contribution. It is roughly consistent with the hole Fermi surface pockets about 
zero 2D momentum revealed by ARPES in the iron-pnictide superconductor BaFe2_ x Co 3: As2 
jsljs]. (See fig. [21 caption.) Figures [2b and [2b show a mix of electron and hole structure 
at cSDW momentum (7r/a)x(y). The "V" shape that separates pink from purple in fig. 
[2b and that separates purple from black in fig. [2b defines an electron Fermi velocity that 
coincides with the cSDW spinwave velocity along the corresponding principal axis: v\ = vq 
and Vt = vq/'-Jcsdw- ARPES measurements reported in ref. sj also reveal electron Fermi 
surface pockets centered at cSDW momenta, with Fermi velocities that show the same 
type of anisotropy. The electron Fermi velocities extracted from the inner edge of the 
"V" in their dispersion curves, in particular, are v\ ~ 0.7 A-eV and v t ~ 0.3 A-eV. These 
values are remarkably close to the corresponding spin-wave velocities extracted from inelastic 
neutron scattering measurements on the same material [3]! Figures [2b and [2b also predict 
stronger hole dispersion at ou — — £/[k — (7r/a)x(y)], however. It is due to the divergence 
of the coherence factor sinh 2 #q in Eq. ([9]) as the spin- wave frequency Ub(q) vanishes at 
q = (7r/a)x(y). Such peaks in ImG(fc,u) must therefore disappear just off the QCP, at 
cSDW gaps A cSDW > Ak peak ■ v . Figures Eb and [2b yield a peak width for the hole 
dispersion of /S.k pea ka ~ 0.1. The hole-dispersion peaks about cSDW momenta are therefore 
fragile, and they may not survive effects such as hole (slave fermion) damping and spin- 
wave (Schwinger boson) damping that are not accounted for within the present mean-field 
approximation. 

We have confirmed the main results of the above mean field theory analysis of the hidden 
half-metal phase by obtaining the exact low-energy spectrum of one hole in the two-orbital 
t-J model (CQ) over a periodic 4x4 lattice of iron atoms using the Lanczos technique 24j. The 
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Hilbert space was confined to the S z = 1/2 subspace, and its dimension was reduced further 
by exploiting translation and reflection symmetries. The Heisenberg exchange terms of the 
Hamiltonian §I§ were stored in memory as permutations of the spin backgrounds. Hole 
hopping was computed at each application of the Hamiltonian (CQ), on the other hand. Last, 
the Lanczos procedure was applied numerically using the ARPACK subroutine library [25)]. 
Figures [3^, [lb and [3b show how the magnetic order and how the low-energy spectrum of 
the t-J model ([1]) evolve with the strength of the Hund's coupling. Figure [3^ displays a 
quantitative match between the dispersion of the exact lowest-energy spin-1/2 excitations 
and the prediction for spin-wave excitations about order at J = 0. Next, fig. 

[Tb shows the exact spectrum at the QCP, where degenerate groundstates exist at momenta 
zero and (7r/a)x(y). The lowest-energy excitations again have spin 1/2, and their dispersion 
is qualitatively similar to the spin-wave prediction displayed by fig. [1^. This quantum- 
critical state shows weak cSDW order and moderate spin order. Last, the Hund's 
coupling used in fig. [lb is the critical value — Jo c = 2.27Jf L . It is considerably larger than 
the corresponding value of —Jq c = 1.35JT 1 obtained in the absence of a hole by one of the 
authors This is consistent with the previous mean field theory result for Jq c , and it is 
likely due to the suppression of quantum fluctuations in the antiferromagnetic state by intra- 
orbital hole motion. These matches indicate that the above Schwinger-boson-slave-fermion 
mean field theory is a valid approximation of the 2-orbital t-J model, Eq. (JTJ, in the hidden 
half-metal phase. Finally, fig. [3b demonstrates that the groundstate is a robust cSD W state 
with a large ordered moment if Hund's rule is obeyed. It carries momentum ±(7r/2a)x(y) 
(cf. ref. 15j), however, and it therefore is unable to account for any of the Fermi surface 
pockets that are observed by ARPES on iron-pnictide materials js ] {7] 8]. 



III. CONCLUSIONS 



Above, we have shown that a hidden half-metal state near a QCP into a cSDW 
state exhibits the nested 2D Fermi surfaces that are characteristic of iron-pnictide high- 
T c superconductors jsj. The Fermi surfaces predicted here with the bare minimum of 3d xz 
and 3gL, orbitals are in fact similar to those obtained by electronic band structure calcula- 

n 

tions that include all five 3d orbitals[5[. In particular, zone-folded Fermi surfaces centered 
at momentum (it /a, tt/o) have low spectral weight in such case [cf. fig. [2H]. Further, our 
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mean field theory predicts a mixture of electron and of hole dispersion for Fermi surface 
pockets centered at cSDW momenta (see fig. [2]). This has been been observed by ARPES 
on iron-pnictide superconductors in certain cases 8] . 

We have also demonstrated above that the low-energy spectrum of the hidden half-metal 
state at the QCP contains zero-energy spin-wave excitations that disperse anisotropically 
away from cSDW momenta. The critical spin-wave spectrum notably exhibits a local maxi- 
mum at the Neel momentum (ir/a, 7r/a), which is consistent with INS on a parent compound 
to iron-pnictide superconductors^. It was pred.cted previously by one of the authors iu 
the case where mobile holes are absent [14fl. Last, we have identified the anisotropic spin- 
wave velocities at cSDW momenta with the corresponding electron Fermi velocities around 
Fermi surface pockets centered at cSDW momenta. A comparison of independent ARPES 
and INS on the same iron-pnictide compound bears out this identification |8| 18]. 

J.P.R. thanks Edward Rezayi, Vitor Vieira, and Maria Jose Calderon for discussions, and 
he thanks Veronique Brouet for correspondence. Exact diagonalization of the t-J model (CQ) 
was carried out on the SGI Altix 4700 at the AFRL DoD Supercomputer Resource Center. 
This work was supported in part by the US Air Force Office of Scientific Research under 
grant no. FA9550-09- 1-0660 and by the FCT under grant PTDC/FIS/101126/2008. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008) 
[2] C. de la Cruz, Q. Huang, J.W. Lynn, J. Li, W. Ratcliff, J.L. Zarestky, H.A. Mook, G.F. Chen, 

J.L. Luo, N.L. Wang and P. Dai, Nature 453, 899 (2008). 
[3] K. Haule, J.H. Shim, and G. Kotliar, Phys. Rev. Lett. 100, 226402 (2008). 
[4] F. Ma, Z.-Y. Lu and T. Xiang, Phys. Rev. B 78, 224517 (2008). 

[5] S. Graser, T.A. Maier, P.J. Hirschfeld and D.J. Scalapino, New J. Phys. 11, 025016 (2009). 
[6] J. Fink, S. Thirupathaiah, R. Ovsyannikov, H. A. Durr, R. Follath, Y. Huang, S. deJong, M. 

S. Golden, Yu-Zhong Zhang, H. O. Jeschke, R. Valenti, C. Felser, S. Dastjani Farahani, M. 

Rotter, and D. Johrendt, Phys. Rev. B 79, 155118 (2009). 
[7] V. Brouet, M. Marsi, B. Mansart, A. Nicolaou, A. Taleb-Ibrahimi, P. Le Fevre, F. Bertran, 

F. Rullier-Albenque, A. Forget, and D. Colson, Phys. Rev. B 80, 165115 (2009). 
[8] V. Brouet, M. Fuglsang Jensen, A. Nicolaou, A. Taleb-Ibrahimi, P. Le Fevre, F. Bertran, A. 



10 



Forget, and D. Colson, larXiv: 1105.56041 . 
[9] J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. Chen, X. Dai, J. L. Luo, 

Z. Fang, N. L. Wang, Euro. Phys. Lett. 83, 27006 (2008). 
[10] Q. Si and E. Abrahams, Phys. Rev. Lett. 101, 076401 (2008). 
[11] J. P. Rodriguez and E.H. Rezayi, Phys. Rev. Lett. 103, 097204 (2009). 
[12] R. Applegate, J. Oitmaa, R. R. P. Singh, Phys. Rev. B 81, 024505 (2010). 
[13] B. Schmidt, M. Siahatgar, P. Thalmeier, Phys. Rev. B 81, 165101 (2010). 
[14] J.P Rodriguez, Phys. Rev. B 82, 014505 (2010). 
[15] C.L. Kane, PA. Lee and N. Read, Phys. Rev. B 39, 6880 (1989). 

[16] Y. Ran, F. Wang, H. Zhai, A. Vishwanath, and D.-H. Lee, Phys. Rev. B 79, 014505 (2009). 
[17] Jun Zhao, D. T. Adroja, Dao-Xin Yao, R. Bewley, Shiliang Li, X. F. Wang, G. Wu, X. H. 

Chen, Jiangping Hu, Pengcheng Dai, Nature Physics 5, 555 (2009). 
[18] C. Lester, Jiun-Haw Chu, J. G. Analytis, T. G. Perring, I. R. Fisher, S.M. Hayden, Phys. 

Rev. B 81, 064505 (2010). 
[19] C. Edmiston and K. Ruedenberg, Rev. Mod. Phys. 35, 457 (1963). 
[20] PW. Anderson, Phys. Rev. 79, 350 (1950). 
[21] A. Auerbach and B. E. Larson, Phys. Rev. B 43, 7800 (1991). 
[22] D.P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 (1988). 
[23] A. Auerbach and D.P. Arovas, Phys. Rev. Lett. 61, 617 (1988). 
[24] C. Lanczos, J. Research National Bureau of Standards 45 (4), 255 (1950). 
[25] R.B. Lehoucq, D.C. Sorensen and C. Yang, ARPACK Users' Guide (SIAM, Philadelphia, 

1998). 



11 



(A) SPINWAVES AT QCP (B) 4x4x2 SPINS + 1 HOLE AT QCP 
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FIG. 1: The mean-field result for the dynamical spin response function in the zero-temperature- 
large-so limit, Eq. ([5]), is evaluated (a) with the following set of parameters: jf = 0, Jf 1 > 0, 
jj = 0.3 Ji = J^~, t\ = —hJi, = 0, x = 0.01, so = 1/2, and Jo = Joe- Low-energy contours are 
displayed in the inset. Also shown (b) is the low-energy spectrum of the corresponding t-J model, 
Eq. (H|), over a 4 x 4 x 2 lattice with one hole at Jo = — 2.27 Jj 1 . Ordered magnetic moments 
over the groundstate (at zero 2D momentum) are also listed there: (n(k) ■ n(—k))o, where n(k) = 
[2fi B /(N Fe - |)]Ei=oEie i(fcoa+k ' ri) Si ia , and where k = (vr,0,0) and (0,vr/a,0), respectively, for 
hidden ferromagnetic (hFM) and for cSDW order. In general, /ip e = (33/31) (2^b) 2 - 
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(A) IMAGINARY G (k,co) (1/j/) at QCP: k = 
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(B) IMAGINARY G (k,a>) (1/j/) at QCP: k = 
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(D) FERMI SURFACE AT QCP 
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FIG. 2: Shown is the imaginary part of expression ([9]) with the parameter set that is listed in 
the previous caption, and in the following limits listed in order: T — > 0, then sq — > oo, and then 
AcSDW 0. The Fe-Fe distance a = 2.8 A yields the Fermi momentum kp = 0.13 A -1 . Setting 



14]) yields a spin- wave velocity Vq ~ 0.6 A-eV, with anisotropy 



sqJi ~ 70 meV (see text and ref. 
parameter jcSDW = 2.64, and a Fermi velocity vp ~ 1.4 A-eV. Also shown are the Fermi surfaces 
predicted by Eq. Q. Each is doubly degenerate because of the 3d xz and 2>d yz orbitals. 
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(A) HIDDEN HALF METAL, NO HUND'S RULE 
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(B) cSDW METAL, STRONG HUND'S RULE 
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FIG. 3: Shown are low-energy spectra of the t-J model, Eq. ([I]), over a 4 x 4 x 2 lattice with one 
hole, with the parameters listed in the text. Each point is doubly degenerate because t± = 0. A 
comparison with the hole spectrum, £f(k), and with the spin- wave spectrum, u>b(k), is made (a) at 
x = in the absence of Hund's rule, Jq = 0. The latter is enforced (b) by setting Jq = — 23 J±. 
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